An Analytical Approach to the Spin Distribution of Dark Halos 
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ABSTRACT 

We derive the distribution of the dimensionless specific angular momentum 
of dark matter halos, P(j), in the framework of the standard tidal torque the- 
ory, and explain the characteristic shape of P(j) commonly observed in N-body 
simulations. A scalar quantity (shear scalar), r, is introduced for measuring the 
effective strength of the tidal torque force acting upon the halo from the sur- 
rounding matter. It is found that the ubiquitous and broad shape of P(j) can 
mostly be attributed to the unique property of the shear scalar r, and also that 
P(j) is insensitve to the underlying collapse dynamics. Our result demonstrates 
that although the shape of P(j) is ubiquitous and close to log- normal, the dis- 
tribution is not exactly log-normal but decays exponentially at the high angular 
momentum end, and drops slowly as a power-law with an index 2 at the low 
angular momentum end. 

Subject headings: cosmology: theory — large-scale structure of universe 

1. INTRODUCTION 

N-body simulations have shown that the dark matter halos of the universe have 
the intrinsic angular momentum, the distribution of which has a broad and ubiquitous 
shape independent of other halo properties such as the halo mass scale, the halo 
formation epoch, and so on (Barnes & Efstathiou 1987; Heavens & Peacock 1988; 
Lemson & Kauffmann 1999; Bullock et al. 2001; Chen & Jing 2002). This characteristic 
distribution of the halo angular momentum is often approximated by the log-normal 
distribution with ubiquitous mean and variance in the log. Although it has been proved 
that this characteristic distribution can be produced by the current standard model 
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(Maller et al. 2001), the following important questions have yet to be answered: what the 
physical origin of the characteristic distribution is, and whether it is truly log-normal. 

The standard model for the origin and evolution of the halo angular momentum is 
the tidal torque theory (Doroshkevich 1970; White 1984). In this model, the angular 
momentum of a proto-halo, L, is generated by the misalignment of the principal axes of the 
local shear tensor with that of the inertia tensor and grows with time in the Lagrangian 
coordinate: Li(t) = —a 2 (i)D(t)€ijkTjiIi k . Here, a(t) is the scaling factor, D(t) the growing 
rate of density perturbations, and hk{— J PoqiQkd 3 c[ with po being the mean density) is the 
inertia tensor of a proto-halo site in the Lagrangian q space, representing the geometrical 
shape of a proto-halo. In addition, T^(= djdi<f> with (f> being the gravitational potential) is 
the gravitational shear tensor that quantifies the tidal torque from the surrounding matter. 

Strictly speaking, the tidal torque theory is a linear theory, and should break 
down in the non-linear regime since it does not account for the complicated nonlinear 
process of halo formation after recollapse. To one's surprise, however, it has been 
demonstrated by N-body simulations that the linear tidal torque theory works well even 
in nonlinear regime (Frenk et al. 1988; Heavens & Peacock 1988; Catelan & Theuns 1996; 
Sugerman et al. 2000; Mailer et al. 2001). Therefore, we will adopt the tidal torque model, 
and derive the distribution of the halo angular momentum. 

2. DERIVATION 

2.1. SEPARATION OF THE VARIABLES 

For the distribution of halo angular momentum, it is more convenient to use a 
dimensionless quantity than to use L itself since L scales as M 5//3 : L ~ I ~ J q 2 d 3 q ~ 
R 5 ~ M 5 / 3 , with R and M being the typical size and mass of a proto-halo respectively. 
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We thus define a time-independent specific angular momentum, j, as ji = e^Tjilik with 
I = IM~ 5 / 3 ), since the log(j)-distribution is evaluated at the same epoch and has a same 
shape as the log(L)-distribution, except for a constant horizontal shift. 

Rotating the coordinates into the local shear principal axis frame, we find that the 
magnitude of the specific angular momentum is given as 

j = {/ 2 3 (A 2 - A 3 ) 2 + J* (Ai - A3) 2 + - \ 2 ) 2 } 1/2 , (1) 

where Ai, A2, A3 are the three eigenvalues of the local shear tensor, and I23, I31, In are the 
three off-diagonal components of the dimensionless inertia tensor, I. It is worth noting 
that only the off-diagonal components of the inertia tensor in the shear principal axes are 
involved in the generation of the halo angular momentum. The diagonal components of 
the inertia tensor are involved in the gravitational collapse. Interestingly enough, however, 
it has been recently found by N-body simulations that the shear and inertia tensors are 
strongly, but not perfectly, aligned (Lee & Pen 2000; Portiani et al. 2001b). In other words, 
in the shear principal axes, the magnitudes of the inertia off-diagonal components are 
small relative to the diagonal parts. This numerical evidence for the small off-diagonal 
components of the inertia tensor prompts us to regard the misalignment to be statistical 
fluctuation and 123,131, In as independent random variables. Now that the three eigevalues 
of the shear tensor are also random variables, and equation (1) describes how the positive 
random variable, j, depends on the random variabes, I and Aj's. 

To derive the distribution of j, we first separate the variables as follows. We 
view the generation of the halo angular momentum as a 3 x 3 diagonal tidal tensor 
Rij = Diag[ri, r i, acting on an inertia-moment vector I, i.e., ji = VSRijlj, in the 
frame of local shear-principal-axis, where the components of and I are given as 
r\ = (A 2 - A 3 )/\/3, r 2 = (A 3 - Ai)/\/3 and r 3 = (Ai - A 2 )/v / 3, and h = I 23 , h = hi 
and I 2 = In- As the three components of I are independent random variables, we may 



- 5 - 



further express them in spherical-polar coordinates as ii = 7 sin cos 0, 7 2 = I cos9, 
7 3 = /sin 9 sin with 7 = (7 2 + 7f + I 2 ) 1 / 2 , 9 = cos^ih/I), = tm-^h/h). On the other 
hand the diagonal components of are not all free but constrained by the condition of 
r± + r 2 + r 3 = 0. One may conveniently regard the diagonal components (ri,r 2 ,r 3 ) also as 
a vector r; this vector lies in a plane perpendicular to the (1, 1, 1) symmetric axis in the A 
space, or r = A x (1, 1, l)/^/3. Therefore, r depends only on one polar angle, say ip. To 
find the polar-coordinate expression for r, we perform a coordinate transformation from 
Ai, A 2 , A 3 into x, y, 5, with x = (Ai - 2A 2 + A 3 )/-\/6, y = (Ai - X 3 )/V2, 5 = A x + A 2 + A 3 . 
Under this transformation, r has only two components: r = (x,y,0). So, we can express 
x = rcos-0, y = rsinifj. Now, rotating the axes back into the shear principal axes, we have 
a polar-angle expression for r in the shear principal axis frame: r\ = cosip + smip, 

r 2 = — ^rsin-0, r 3 = cosip + sin-0. 

Using these angular variables for I and r, we find the following expression for j with 
separated variables: 

j = IrY^,9A), (2) 

where 

r 1 1 1 2 -] 1 / 2 

Y(tjj, 9, 0) = < sin 2 9(- cos 2 ip + - sin 2 if) -= cos 20 sin 2ip) + - cos 2 9 sin 2 1[) > . (3) 

1 2 6 2y3 3 J 

The off-diagonal components of the inertia tensor Ii, J 2 and J 3 are further assumed to be 

Gaussian variables (Catelan & Theuns 1996). So, the distribution of / = \J 1 2 + If + if is 

a weighted Gaussian given as = v^a7 3 / 2 exp (-/ 2 /2 < t|) (if, where 07 is the 

standard deviation of each off-diagonal component. While the solid angles, and 9 are 

uniformly distributed on the sphere. 

Though the shear scalar r was first pointed out to play an important role for non- 
spherical halo formation and angular momentum generation by Chiueh & Lee (2001), Sheth 
& Tormen (2002) recognized that r can in fact be conveniently expressed as [^XlLi^] 1 ^ 2 
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where y±, ■ ■ ■ ,y$ are mutually independent Gaussian variables, having the same standard 
deviations equal to that of the density field, as- Therefore, the distribution of r is also a 
weighted Gaussian: P r (r) = v /o7^(25/12cr|)r 4 exp (-5r 2 /4cr|). Sheth & Tormen (2002) 
commented that the broadness of P(j) may be attributed to the shape of P r (r). However, 
P(j) is contributed by the conditional distribution of r only for those regions out of which 
dark halos condense, i.e., P r (r | halo), instead of the unconditional P r (r) given above. In 
the Press-Schechter model of halo formation (Press & Schechter 1974) where the collapse 
threshold depends only on 5 and not on r, P r (r |halo) indeed equals P(r). However, the 
actual halo collapse is generically non-spherical, which leads to higher collapse thresholds 
depending both on 5 and r (Sheth et al. 2001; Chiueh & Lee 2001). This fact differentiates 
the actual halo r-distribution per mass bin, P r (r |halo), from the unconditional P r (r). 

To find P r (r | halo), we examine the non-spherical collapse using the collapse condition 
given by Chiueh & Lee (2001): £ = (l + ^ ■ Chiueh & Lee (2001) originally suggested 
a — P — 0.15; later Lin et al. (2001) refined the values as a = 0.26 and (3 = 0.16, which are 
adopted here. Since the collapse condition depends on r, P r (r|halo) should no longer equal 
P r {r). We calculated P r (r|halo) numerically using the random- walk method described in 
Lin et al. (2001). Two non-trivial results are found. First, P r (r|halo) is independent of mass 
up to the numerical accuracy, and second, P r (r) turns out to still be a good approximation 
to P r (r | halo). Figure 1. compares the numerical P r (r|halo) near the characteristic mass 
M, with P r (r). The good argreement suggests that the distribution of j is not sensitive to 
the detailed dynamics of halo formation. This feature is the origin of the scale-independent 
shape of the halo angular momentum distribution to be shown below. 
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2.2. The j-Distribution 

Now, replace P r (r|halo) by the simpler P r (r). By equation (2), the angular distribution 
P(j) can be derived as follows: 

P(j) = J dX^dXs J dIsm9d9d(j>P x (\ 1: \ 2: \ 3 )P i 5 D (j - IrY{9,<f>,^)), (4) 

where 

Px(\i, A 2 , A3) = exp (-^ - ^ J | (Ax - A 2 )(A 2 - A 3 )(A 1 - A 3 )|, (5) 

and the Dirac ^-function constrains j to other variables through eq.(2). Note that the last 
factor in eq.(5) is no more than |rir 2 r 3 |(= J(r,ip)), arising from the Jacobian of the angular 
distribution of the shear principal axes(Doroshkevich 1970). Thus, J(r,ip) = r 3 J(ip), with 
J(ip) = | (3 cos 2 ip — sin 2 -?/') sin tp\ . The angular dependence of J(ip) results in P\ to contain 
three null lines, separated by 60 degrees, on the plane where the r vector lies. The null lines 
correspond to the degeneracies, Ai = A 2 , A 2 = A3 and A3 = Ai. This feature renders the 
random angle ip to distribute non-uniformly with three-fold symmetry. 

Upon changing the variables from Ai, A 2 , A 3 to 5,r,ip, the 5, I and r integrals can all 
be evaluated analytically, and eq.(4) now becomes 

PV)=*>- 1 3 3 f ^ K ^W )m (6) 

where d 3 Q = sin OdOdifidip, K x is the lst-order modified Bessel function, b = y/ /asoJ, w the 
normalization factor and Y is given in eq.(3). As log(j'/6) = log(j') — log(fe), the shape of 
the log(j)-distribution is not affected by the value of the unknown b although the peak 
position is. The angular integration of eq.(6) can only be evaluated numerically and the 
solid line of Figure 2 presents the resulting jP(j). It shows that though the main body of 
P(j) resembles the log-normal distribution (dashed line), the distribution is not exactly 
log-normal. One may check from eq.(6) by asymptotic expansion that P(j) drops off 
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exponentially at the high angular momentum limit (cf. Mailer et al. 2001) and behaves 
as a power law (~ j 2 ) in the lower angular momentum limit. The best-fit log-normal to 
our distribution has a log- width ciogj = 0.63 (dashed line), which is in fair agreement 
with the log-width, <7i og j = 0.5, obtained from N-body simulations reported by Bullock et 
al. (Bullock et al. 2001). (Their definition of A' is equivalent to our j modulo a constant 
factor.) 

For comparison, we have also simulated non-spherical halo collapse using the random- 
walk algorithm. When the collapse condition given in the previous section is satisfied, the 
shear tensor is diagonalized (Chiueh and Lee, 2001) and the three principal eigenvalues \'s 
are used to evaluate the halo shear vector r. The angular momentum of the dark halo is 
then computed according to eq.(l). Figure 2 also plots the resulting jP(j) (solid squares). 
To ensure the accuracy of the our method, we have also repeated the same procedure for 
the Press-Schechter model and computed its jP(j) (solid triangles), to which the solid 
line should be identical. As one can see, all three results agree with each other, revealing 
that the angular momentum distribution is insensitive to the detailed dynamics of halo 
formation. 

We would like to stress that the unique shape and log(j') width of the angular- 
momentum distribution are the direct outcome of our model. It involves no fitting 
parameter. 

3. DISCUSSIONS AND CONCLUSIONS 

To explain the observed characteristic shape for the angular momentum distribution of 
dark halos in N-body simulations, we have adopted the standard tidal torque theory to study 
the proto-halo angular momentum j iy which can succinctly be represented by a random 
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shear acting upon a random vector Jj. The diagonalized Rij, whose three components 
are the mutual differences of the three eigenvalues of the local shear tensor, quantifies 
the non-sphericality of the local gravitational potential, and the vector I characterizes 
the misalignment between the matter distribution and local gravitational potential. With 
this picture, we have derived a new expression for the halo angular momentum, which 
separates the variables into three parts: r, J, Y(ip, 9, <fi). With P r (r|halo) being shown to be 
approximately a weighted Gaussian and assuming Pj(I) to also be a weighted Gaussian, we 
have derived the angular-momentum distribution, i.e., P(j), to be close to log-normal, with 
a mass-independent width in fair agreement with that determined from simulations. This 
"quasi" -log- normal distribution is simply a consequence of the nonlinear coupling of five 
random variables, r, /, 9, and ip. 

Though the distribution of log(j) has a mass-independent shape, the mean of 
log(j) is predicted to depend on b (= aj(Js(M)). As our framework does not permit 
the mass dependence of 07 to be determined, this present work can only account 
for the ubiquitous shape, but not the peak location, of the angular momentum 
distribution. Nonetheless, simulations have empirically found that even the mean value 
of log j is also mass independent (Barnes & Efstathiou 1987; Heavens & Peacock 1988; 
Lemson & Kauffmann 1999; Bullock et al. 2001; Chen & Jing 2002). Incorporating this 
finding, we thus demand that 

aj{M)^Gj\M). (7) 

Such mass dependence of aj is turned to be the prediction of the present work. That is, the 
misalignment (per mass) I has a wider distribution for larger halos than for smaller halos at 
the same epoch, since as is greater for smaller halos than for larger halos. This prediction 
does make qualitative sense, in that at a given epoch large halos have undergone vigorous 
growth with mergers of sub-halos; it thus makes the interior matter distribution of large 
halos less correlated with the surrounding matter from which the gravitational forces are 



-10- 



exerted. By contrast, small halos have grown mildly and hence tend to be in a quiescent 
state, thus keeping themselves well correlated with the surrounding matter. At any rate, 
this prediction as well as the Gaussianity assumption of / are both testable by cosmological 
simulations. 

This work has been supported by the Taida-ASIAA CosPA Project. T.C. acknowledges 
the partial support from the National Science Council of Taiwan under the grant: 
NSC90-2112-M-002-026. 
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Fig. 1. — The distribution of the halo shear scalar r in the logarithmic scale. The solid 
squares are P r (r|halo) near M* obtained from the Monte-Carlo simulation using the non- 
spherical collapse condition. The solid curve is P r (r) given by the analytic formula with the 
Press-Schechter collapse condition. 
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Fig. 2. — The distribution of the logarithmic specific angular momentum of dark halos, 
jP(j). The solid curve is our analytic result given by eq.(6); the squares are Monte Carlo 
result for the non-spherical dynamical model, and the triangles are Monte Carlo result for 
the Press-Schechter model. The dashed line is the best log-normal fit to jP(j). 



